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1 Introduction 



The recent discovery of anisotropy in the cosmic microwave background (CMB) lends support 
to the hypothesis that structure in the universe formed via gravitational instability from small 
density perturbations. Furthermore, the observed anisotropy allows us to characterize the initial 
perturbations statistically with enough precision to draw important cosmological conclusions. In 
particular, the COBE detection (Smoot et al. 1992, Bennett et al. 1994) shows that density per- 
turbations existed on scales larger than the horizon at the epoch of recombination. Furthermore, 
the data are consistent with the hypothesis that the perturbations are Gaussian distributed and 
can be used to place constraints on the power spectrum of the initial density perturbations. CMB 
anisotropy is already one of the most important pieces of cosmological data, and in coming years 
its importance will only increase. 

In recent years, many cosmologists have favored theories in which the density parameter f^o is 
equal to one. In particular, this value for f^o is predicted by the simplest versions of the popular 
inflationary scenario. However, many dynamical measurements indicate a low-density universe 
(Peebles 1993, Ratra & Peebles 1994a). It is well known that the standard COBE- normalized 
Cold Dark Matter (CDM) scenario, in which universe is flat and contains only baryons and 
CDM particles, and in which the initial density perturbations are of the Harrison-Zel'dovich 
type, predicts fluctuations that are too large in amplitude on scales of galaxy clusters and 
below, although a number of slight variants on this model can be devised that fit the data better 
(e.g., Bunn, Scott, & White 1995; White et al. 1995). 

There are two classes of low-density cosmological model that can be motivated by the in- 
flationary universe scenario. One is the A-model. In this model the matter density is small, so 
0,q < 1; however, the cosmological constant contributes to the total mean density required to 
make the universe flat: f^o + ^a = 1- The usual inflationary scenarios predict that the universe 
is flat, so that f2o + ^A must be extremely close to unity (Kashlinsky, Tkachev, & Frieman 1994). 
The cosmological constant is equal to the vacuum energy density of the Universe. Although we 
have no reason to be certain that it must be zero, the values that are of interest to cosmologists 
(i.e., the values that make Oa of order unity) are unnaturally small from the point of view of 
particle physics (Weinberg 1989), and so these A- models are often regarded as unappealing. 

The other class of low-density inflationary scenario is the open model, in which the universe 
has negative spatial curvature. Recently, the possibility of realizing an open universe has been 
discussed in the context of inflation theory (Bucher, Goldhaber, & Turok 1994; Yamamoto, 
Sasaki, & Tanaka 1995; Linde 1995; Linde & Mezhlumian 1995). The essential idea is based 
on the semiclassical picture of a bubble nucleation, which is described by a bounce solution 
(Coleman & De Luccia 1980). One bubble nucleation process can be regarded as the creation 
of a homogeneous and isotropic spacetime with negative spatial curvature inside the bubble due 
to the 0(4)-symmetry of the bounce solution. 

It is of great interest to ask whether an open universe created in this one-bubble inflationary 
scenario is observationally acceptable or not. Open CDM models have been investigated by many 
authors (Lyth h Stewart 1990; Ratra & Peebles 1994ab; Sugiyama & Silk 1994; Kamionkowski 
et al. 1994; Gorski et al. 1995; Liddle et al. 1995). Their investigations are based on the simple 
assumption that the quantum state of a scalar field is in the conformal vacuum state at the 
inflationary stage; however, this is unlikely to be the prediction of the one-bubble inflationary 
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scenario. It has been pointed out that if we take the Bunch-Davies vacuum state as the state 
of scalar field, the CMB anisotropy in a low fio universe appears quite different due to the 
super-curvature mode (Yamamoto, Sasaki, & Tanaka 1995). 

Following the usual inflationary picture in which the quantum fluctuation of a scalar field 
generates the density perturbation, we must investigate the quantum state of fields inside the 
bubble. Attempts have been made to study this problem by developing a field-theoretical for- 
malism based on a multidimensional tunneling wave function (Tanaka, Sasaki, & Yamamoto 
1994; Sasaki et al. 1994; Tanaka & Sasaki 1994; Yamamoto, Tanaka, & Sasaki 1995; Hamazaki 
et al. 1995). Bucher et al. have also considered this problem (Bucher, Goldhaber & Turok 1994; 
Bucher & Turok 1995). However, this problem requires further investigation. 

In this paper, we consider the simple case in which the quantum state of a scalar field is 
in the Bunch-Davies vacuum state, and compare the predictions of this model with several 
cosmological observations. This situation is physically definite and clear, and as long as the 
bubble nucleation occurs in the de Sitter inflationary background (where the initial inflationary 
period has lasted sufficiently long), the field can be approximated by the Bunch-Davies vacuum 
state, provided that the effect of the bubble nucleation process is negligible. 

In section 2, we first consider the initial spectrum by solving for the evolution of cosmological 
perturbations in the open inflationary stage. In section 3 we use this initial power spectrum to 
calculate various cosmological quantities and compare these predictions with observations. We 
also compare these results with those of a previous analysis of the open inflationary model based 
on the conformal vacuum state, and a A-model with Harrison-Zel'dovich spectrum. Section 4 is 
devoted to a discussion of our results. We will work in units where c = 1 and Ti = 1. 



2 Initial Conditions 

In this section, we consider the evolution of cosmological perturbations in an open inflationary 
stage and derive the initial spectrum of perturbations. Ratra & Peebles (1994b) have investigated 
the evolution of cosmological perturbations in an open inflationary universe with gauge fixed. 
Bucher, Goldhaber h Turok (1994) have also investigated cosmological perturbations in an 
inflationary stage in the gauge-invariant formalism. In the first half of this section, we follow 
the work by Bucher, Goldhaber & Turok (1994). 

In an open universe the line element can be written as 



ds 2 = a 2 (r?) 



—drf + dx 2 + sinh 2 x<^(2) 



(1) 



where rj is conformal time. We consider a scalar field <p with potential V ((/>), which drives 
inflation. Writing the scalar field as a homogeneous part and a small inhomogeneous part, i.e., 
(j) = 4>q + 5(fr, the equations for the homogeneous part are 



bl + 2H<t>' + a 2 ^ = 



and 



8ttG 



-& + a 2 V 



(2) 
(3) 
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where Ti := a' /a and the prime denotes differentiation with respect to rj. 

The small fluctuation 6(f) gives rise to a metric perturbation. As we are interested in a scalar 
perturbation in a scalar field dominated universe, the metric perturbation can be written, 



ds 2 = a 2 (rj) 



-(1 + 2<S>)dr] 2 + (1 - 2$) ( d X 2 + sinh 2 X dfl 2 



(2) 



(4) 



corresponds to the curvature perturbation or gravitational potential. Then we get the evolu- 
tion equations for the perturbation (Mukhanov, Feldman, &: Brandenberger 1992), 

+ 2 f H - fy) $' + f-L 2 + 4 + 2H' - 2H^-) $ = 0, (5) 
' + (6) 



where 

and is the Laplacian on the unit sphere. 

Now let us consider an inflationary stage of the universe inside a bubble. To solve the above 
equations analytically, we assume that the potential is nearly flat, and use the approximation, 
V — Vo+V'4>, where V = const. We also assume that the background spacetime is approximated 
by de Sitter spacetime, that is, a(t) = —l/Hsinhrj. 

Then the field equation for the homogeneous part is 



-V 

H 2 sinh 2 T] ' 



% - 2coth #' = tt0 - ]2 , (8) 



with H 2 = 8ttGV /3. 

Eq. (JsT) can be integrated, giving 



. — V' — cosh rj + 3 cosh rj sinh r/ + 2sinh rj 



H 2 3 sinh rj 

and the perturbation equation (||) reduces to, 



(9) 



*"-^*'+H + -f^)*^ <-) 

To solve these equations, we need the initial values. To determine the initial values, we must 
investigate the problem of what the quantum state is inside the bubble. As mentioned before, 
this is a very important problem, which demands further investigation. We here consider the 
case in which the scalar field is in the Bunch-Davies vacuum state. This is the case provided 
that the effect of bubble nucleation process is small and negligible. But we should keep in mind 
that this point needs examination in the various models of one-bubble inflation scenario, taking 
into consideration the effect of bubble nucleation. 
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Recently, quantum field theory in de Sitter space-time associated with the open chart has 
been investigated (Sasaki, Tanaka, & Yamamoto 1995). According to this analysis, a quantized 
scalar field with mass m 2 <C H 2 in the Bunch-Davies Vacuum state is described in the second 
quantized manner as, 



+ h.c, (11) 



(T,Lm 



Irn 



where 



-1 



y/8p(p 2 + 1) sinhvrp 

e np/2 {ip + coth^e^ + ae-^ p/2 {ip - cothr})e ipri 



:i2) 



Pi L (cosh x) 



H 



w WI (x)iwn ( 3)), 



(J > o), 



(13) 



<r takes on the values ±1, a is the annihilation operator, and T(z) is the gamma function. The 
orthonormal harmonics on a three-dimensional unit hyperboloid Y p i m (x,Q) are 



Yplm(X> ^ (2)] 



pT(ip + 1 + 1) 



T(ip+1) 



P ip-l/2( COsh X) 



\J sinh x 



Yi m [Q.(2))i 



(14) 



with normalization 



dx / dfi (2 )sinh X^imi(x,^(2))^ 2 i2m 2 (x,^(2)) = <*(Pi -^2)^/2^1 



m 2 • 



(15) 



In the above expression, the usual harmonics behave as Y p i m oc e~ x at scales larger than the 
curvature scale, x>l, although vu\i m is constant for x 3* 1. Thus v^y m represents a fluctua- 
tion larger than the curvature scale, so we call this mode a super-curvature mode. The necessity 
of super-curvature modes for a complete description of a random field in an open universe has 
also been discussed by Lyth & Woszczyna (1995). 

Next, let us consider the curvature perturbation, which can be written in the mode-expanded 
form 

roc 

$ = y dp ]T $ Piff ( J 7)y pim (x ) n ( 2)) + E*(*)^) w (*)iW y im(fi(2))- (16) 

^ a,l,m l,m 

For the continuous mode (p,l,m), Eq.(||) reduces to 



p 



6(1 



3-e 27 ? 



4(3 + e 2 ^) 
3 - e 2? ? 



$ p = 0. 



(17) 
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The solution that behaves like & p — > e^ 1 * P ' JJ as r\ — > — oo is (Bucher & Turok 1995) 

•>-"-*(i+£5t)- < 18 > 

The equation for the super-curvature mode is 

6(1 -e^) f 4(3+g>) \ _ 

*(*)" 3 _ e 2r ? $ (*) + ^ 4 ~ 3 _ e 2^ J $ (*)-°» ( 19 ) 

and we find a solution 

£(,) = e 2 ". (20) 

To determine the amplitude of we use Eq.(^). From the behavior at r/ — > — oo, we find 
for the continuous mode 

Wttf = ; 1 L^^$ p + ae-^±±^*X (21) 

p ' K " H V8p(P 2 + l)sinhvrpl 2 - ip p T 2 + ip P J' K J 

and for the super-curvature mode 

^)(l)~^ W (^ (22) 
We therefore have the following spectrum at the end of inflation, by taking the limit 77 — ^ 0, 
d> p (0) 2 := lim *pAv? 



a=±l 

/2ttGV'\ 2 cothvrp p 2 + 1 
V -H" / 2p(p 2 + 1) p 2 + 4 j?™ 
MvrGy'X 2 cothvrp 



and 



\ 3H J 2p(p 2 + l] 



d> w (0) 2 = lim$ w (r ? ) 2 
4ttGV\ 2 1 



(23) 



( •■ > / / y 1 



(24) 



For comparison, we also investigate the case when the scalar field is assumed to be in the 
conformal vacuum state (Lyth & Stewart 1990, Ratra & Peebles 1994b). In this case, the scalar 
field is written 

6(j) = j dpJ2xp{v) Y pim{x,^(2))b p im + h.c. , (25) 

" l,m 

where 

(ip + cothr])e- ip7 i 1 

x " < " ) = —jmrrrwr (26> 
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After a similar analysis, we get the spectrum of curvature perturbations at the end of inflation, 



^W=f^) Twl^-TT- (27) 



tj-o py " \ 3H J 2p(p 2 + l)' 

The conformal vacuum case differs from the Bunch-Davies vacuum case in two ways, the factor 
coth7rp and the super-curvature mode. 

Lyth & Stewart (1990) have investigated perturbations in an open inflationary universe, 
and have given a relation to relate the curvature perturbation and the scalar field perturbation 
TZ ~ —(H/(j))8(j), though a paper justifying this relation has never been published. But the above 
investigation shows the correctness of their result on all scales except for the small difference of 
the former coefficient. 



3 Observational Confrontations 

Now we start testing the predictions of the open universe in the context of CDM cosmology 
with the initial conditions obtained above. The matter-dominated open universe has the line 
element ([[]) with a{rf) = cosh 77 — 1. In this section we use r/(> 0) as the conformal time in the 
matter-dominated universe. 

(1) CMB Anisotropics 

Let us first consider the CMB temperature fluctuation. Having obtained an initial perturba- 
tion spectrum, we can compute the temperature fluctuations in the gauge-invariant formalism 
(Sugiyama & Gouda 1992). As usual, we write the temperature autocorrelation in the form, 

= ^-y2(2l + l)CiPi(cosa). (28) 

Figure 1(a) shows the power spectrum of temperature fluctuations, 1(1 + 1)C; x 10 10 /27r, for 
various values of Qq with initial conditions associated with the Bunch-Davies vacuum state. We 
have taken n B h 2 = 0.0125 and Hubble parameter h = 0.75, 0.70, 0.65, 0.65, 0.60, for Q = 0.1, 
0.2, 0.3, 0.4, 0.5, respectively, to take the age problem into consideration. Figure 1(b) shows 
the same quantities with the conformal vacuum state (Kamionkowski et al. 1994; Gorski et 
al. 1995). For reference, we show the corresponding results for a A-model with a Harrison 
Zel'dovich spectrum in Figure 1(c) (Sugiyama 1995). Here the parameter fish 2 = 0.0128 and 
h = 0.8. We also show the results of several CMB experiments, taken from the paper by Scott, 
Silk & White (1995). Open models may have trouble fitting the data near the "Doppler peak" 
on degree scales, although assessing the significance of this problem will require very careful 
investigation (Ratra et al. 1995). 

The differences between Fig. 1(a) and Fig. 1(b) at low multipoles come almost entirely from 
the contribution of the super-curvature mode (Yamamoto, Sasaki, & Tanaka 1995), 

C w , = {^^-) 2 {\f{VLs)W {if)l ( m - r, LS ) dri'^p-Wwivo ~ (29) 
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where 77^5 and r/o are the recombination time and the present time, respectively, Wuv is defined 
in Eg . (|l3|) , and f(rj) is the decay factor of the curvature perturbation, 

sinh 2 77 — 3ri sinh rj + 4 cosh rj — 4 
/(??) = 5 (cosh ,7 -1)3 • (30) 

The most accurate and reliable CMB anisotropy data at the present time come from the 
COBE DMR experiment. In addition to providing us with accurate estimates of the fluctuation 
amplitude, the data from this experiment can be used to constrain the shape of the power 
spectrum. We have used the two-year COBE data (Bennett et al. 1994) to place constraints 
on open inflationary models, following a procedure based on the Karhunen-Loeve transform 
(Bunn, Scott, & White 1994; Bunn & Sugiyama 1995; White & Bunn 1995; Bunn 1995ab). 
This procedure gives results that are generally consistent with the spherical-harmonic technique 
devised by Gorski (1994). We will now describe this procedure. 

In inflationary cosmological models, the CMB anisotropy is a realization of a Gaussian 
random field. If we expand the anisotropy in spherical harmonics, 

00 / 

AT(r) = J2 E aimY lm (r), (31) 

1=2 m=-l 

then each coefficient a\ m is an independent Gaussian random variable of zero mean. Furthermore, 
the variance of a\ m is simply C\. With this information, we can in principle compute the 
probability density p(d \ C{) of getting the actual COBE data d given a power spectrum Cf. 
since each data point is a linear combination of Gaussian random variables, the probability 
distribution d is simply a multivariate Gaussian, 

p(d\ Ci) oc exp (—d^M^dj , (32) 

where the covariance matrix M can be written in terms of the power spectrum and the noise 
covariance matrix. 

Let us restrict our attention to a few-parameter family of possible power spectra. We will 
denote the parameters generically by q. In this paper, for example, we will consider power spectra 
that are parameterized by two parameters, the density parameter f^o and the power spectrum 
normalization Q = y / 5C^/47r, and so q will be a two-dimensional vector. The probability density 



in equation (32) is then simply the probability density p(d \ q) of the data d given the parameters 
q. If we adopt a Bayesian view of statistics, we can convert this into a probability density for 
the parameters given the data: 

p(q\ d) oc p{d\ q)p(q), (33) 

where p{q) is the prior probability density we choose to adopt. p(q \ d) is generally denoted L(q) 
and called the likelihood. 

The choice of prior distribution is a notoriously troublesome issue. In practice, one generally 
chooses a prior that is a smooth, slowly-varying function of the parameters. In this paper, 
we will adopt a prior that is uniform in fio and one that is uniform in InQ. (This prior is 
approximately equivalent to one that is uniform in the power spectrum normalization C\q near 
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the "pivot point." It differs slightly from one that is uniform in Q, although not enough to affect 
our results significantly.) 

Unfortunately, in order to compute the probability density p(d \ q), and hence the likelihood 
L, one must invert a matrix of dimension equal to the number of data points. For the COBE 
DMR data, this number is of order 4000. Such exact likelihoods have been computed for a small 
class of models (Tegmark & Bunn 1995); however, this is quite a time-consuming procedure. 
The Karhunen-Loeve transform allows us to "compress" the data from 4000 numbers to only 
400 in a way that throws away very little of the actual cosmological signal. The likelihoods 
estimated from the transformed data approximate the true likelihoods well, and are much more 
efficient to compute. For details on how the Karhunen-Loeve transform is performed, see White 
& Bunn (1995) and Bunn (1995ab). 

Once we know L, it is quite easy to place constraints on the parameters q. Since L is a 
probability distribution for q, it should be normalized so that 

jL(q)dq=l. (34) 

Now suppose that we choose some subset R of possible parameter values. Then if 

/ L(q)dq = c (35) 
Jr 

then we can say that q lies in the region R with probability c. If we want to find a 95% confidence 
interval, we simply find a region R such that c = 0.95. One frequently chooses R to be the region 
enclosed by a contour of constant likelihood. 

If one of the parameters is deemed to be uninteresting, the standard practice is to "marginal- 
ize" over it. For example, if we are interested in constraining f2o but not Q, then we replace 
L{Sh,Q) by 

Wg(^o) = J L(n ,Q)dQ. (36) 

This is a natural thing to do: if L is the joint probability density for and then -Z^marg is 
the probability density for alone. 

Figure 2 shows the contours of the likelihood L(flo, Q) for open models associated with the 
Bunch-Davies vacuum state. In computing these likelihoods, we use a linear combination of 
the 53 and 90 GHz two-year COBE maps, with weights chosen to minimize the noise. We use 
the ecliptic-projected maps; maps that were made in Galactic coordinates, and therefore have 
different pixelization, give normalizations that are generally lower by a few percent (Stompor, 
Gorski h Banday 1995; Bunn 1995a). The choice of pixelization appears to affect primarily 
the overall normalization of models; likelihood ratios of models with power spectra of different 
shapes are less affected (Bunn 1995ab). 

Figure 3 shows the marginal likelihoods for Qq for both the Bunch-Davies and the conformal 
vacuum open models. In the Bunch-Davies case, we find that Qq > 0.34 at 95% confidence and 
f2o > 0.15 at 99% confidence. We also show the confidence levels for various Qq hi Table 1. 
For the conformal vacuum models, the likelihood is bimodal, and so the allowed regions are not 
connected. If we take a cut at small f&o and only consider the region > 0.03, we can state 
that at 95% confidence either < 0.085 or fio > 0.36, and 99% confidence, either < 0.14 
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or f^o > 0.23. There are difficulties associated with the interpretation of the likelihoods in this 
case (Gorski et al. 1995). 

Of course, the likelihood L(Qq,Q) provides us with accurate normalizations in addition to 
shape constraints. For any particular value of Oo, we find the value of Q that maximizes 
the likelihood and use this value as the power spectrum normalization. The normalizations 
determined in this way have typical one-sigma fractional uncertainties of approximately 7.5%. 
The maximum-likelihood normalizations computed in this way are listed in the second column 
of Table 2(a) for the Bunch-Davies vacuum case. For comparison, we have also computed the 
conformal vacuum case; these normalizations are given in the second column of Table 2(b). 



(2) Linear density power spectrum 

We next consider the matter inhomogeneities using the COBE normalization as described above. 
As the density perturbation A is related to the curvature perturbation $ by the gravitational 
Poisson equation (Kodama & Sasaki 1984), 

(p 2 + 4)$ p (t/) = 4vrGp(r ? )a 2 (7 ? )A p (r ? ), (37) 

in linear perturbation theory, we can write the power spectrum of the matter perturbation in 
an open universe from Eq.(|23|), 

alP(k) (:= ogA?) = (^^) 2 (p 2 + 4) 2 a^(0) / 2( r?0 ) T( A ; ) 2 

== ^ 2 + 4) 2 ^f y T W 2 , (38) 

where p = a^k, ag = V-^ov 7 ! — ^o> an d Ho = 100/ikm/s/Mpc. The super-curvature mode 
does not contribute on small scales. The model based on the conformal vacuum state leads to 
the same form but without the factor cothvrp. In the CDM cosmology, the following transfer 
function is useful (Bardeen et al. 1986, Sugiyama 1995), 



log(l + 2.34g) 

with 



T{k) = ~ &VJL ' 1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 l 1/4 , (39) 



( — =-] 7= /iMpc^ 1 . 40) 

V 2.7 J Q hexp(-n B - V2h^ B /^o) 



This transfer function is for a flat universe; however, since we are interested in small-scale 
perturbations, the curvature of the universe can be neglected and the transfer function above is 
acceptable. 

The COBE DMR normalization determines the amplitude of the fluctuation. We give the 
numerical value of A determined from the likelihood normalization in the second column in 
Table 2(a). The second column in Table 2(b) gives the value for the conformal vacuum case, in 
which the power spectrum is obtained by Eq.(|38|) without the factor coth7rp. 

We show in Figure 4(a) the density power spectrum a^P^k), for various Qq, with initial 
conditions based on the Bunch-Davies vacuum state. The Hubble parameter is the same as that 
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in Fig.l(a). The points are from Peacock & Dodds (1994). Figure 4(b) shows the density power 
spectra for A-models with Harrison-Zel'dovich power spectrum, in which the Hubble parameter 
is same as that in Fig. 1(c). 

Given the density perturbation spectrum P(k), we are able to calculate erf =((5M/M)g /i _ 1Mpc , 
the variance of the mass fluctuation in a sphere of a radius R = 8/j _1 Mpc, 

a 2 (R) = J k 2 dkP{k)W 2 (kR), (41) 

where the top-hat window function W is defined by W(x) = 3(sinx — xcosx)/x 3 . In Tables 
2(a) and 2(b), we give as for various and h in the open model associated with Bunch-Davies 
vacuum and the conformal vacuum respectively. We also show values for the A model in Table 
2(c) (Sugiyama 1995; Stompor et al. 1995). 

The difference between the Bunch-Davies case and the conformal vacuum case is very small. 
The difference is 10 percent at Qq = 0.05, but is a few percent even at f^o = 0.1. This is because 
the COBE normalization based on a likelihood analysis gives more weight to the behavior of the 
power spectrum at I ~ 10, and less weight at lower multipoles (White & Bunn 1995). When we 
take a <t(10°) normalization, there is a 10 percent difference between the two cases at £lo = 0.1. 
The values obtained are consistent with those in Gorski et al. (1995). 

A precise comparison between predictions and observations of the matter power spectrum is 
difficult. One of the primary problems is that we do not know whether the galaxy distribution 
is an unbiased tracer of the mass distribution. However, if we make the reasonable assumptions 
that the galaxies are not anti-biased and are not extremely strongly biased (say b = a^ 1 Ss 2.5), 
then these calculations suggest that 0.3 ^> Q,q ^> 0.5. Note that the values of f^o preferred by the 
COBE likelihood analysis tend to be higher than this range; however, one might be inclined to 
argue that a model with SIq — 0.4 — 0.5 passes both tests. 



(3)Large-scale bulk velocity 
Next, we consider large-scale bulk velocities, which are given by the following expression, 

4 = / dkP{k)W(kR? exp(-k 2 R 2 ), (42) 

where W{kR) is the window function, and R s = 12/i _1 Mpc is the Gaussian smoothing length 
for comparison with the observational data. In Table 3, we have summarized the computation 
of v 2 R with R = 40/i~ 1 Mpc for various for open models with the initial conditions associated 
with the Bunch-Davies vacuum state. These results are consistent with those of Gorski et al. 
(1995). 

We can compare this results with the recent data from the POTENT analysis (Dekel 1994; 
Liddle et al. 1995): v R=i0h -i Mpc = 373 ± 50km/s. Large values of f^o clearly provide a better fit 
to the velocities. It appears difficult to reconcile models with f^o ^ 0.3 with these data; however, 
it is difficult to make precise statistical statements based on these observations. As is discussed 
by Liddle et al. (1995), this measurement of the bulk velocity contains additional uncertainty 
due to cosmic variance. In addition, it is quite difficult to assess the uncertainties and potential 
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biases in the POTENT analysis, and one should therefore be reluctant to draw firm conclusions 
on the basis of such a comparison. 

(4) epoch of galaxy formation 

Liddle et al. (1995) have performed a detailed investigation of abundances of galaxy clusters 
and damped Lyman-alpha systems in open CDM models using Press-Schechter theory. In this 
paper we will rely on a simple and rough estimate of the epoch of galaxy formation, following the 
work of Gottlober, Miichet, & Starobinsky (1994) and Peter, Polarski, & Starobinsky (1994). 
According to Press-Schechter theory, the fraction of the matter in the universe which is in 
gravitationally bound objects above a given mass Mr at a redshift z has the form 

F(>M « ) = erfc (7I^)' <43> 

where 

= '(%Tl7ta55j- (44) 

where Mr = (4/3)irR 3 p, and f(rj) is the decay factor of the curvature perturbation. 

The choice of 5 C depends on the collapse model. The spherical collapse of a top-hat pertur- 
bation gives S c = 1.69, although non-spherical collapse models suggest other values. Here let 
us consider the range (1.33 < 5 C < 2) (Gottlober et al. 1994). Observations suggest that many 
galaxies seems to have formed at z = 1, then, assuming F{> 1O 12 M ) £ 0.1 at z = 1, we have 
a(M R = 1O 12 M , z = 1) £ 2 ± 0.4. 

Figure 5(a) shows a contour plot of u{Mr = 10 12 M Q , z = 1) in the Qo — h plane for the open 
model. Figure 5(b) is same but for the A- model with a Harrison Zel'dovich spectrum. If we take 
the age problem into consideration, it indicates a lower bound CIq *5 0.4 for the open model. 
Note that the above estimate is very rough, although a similar constraint has been obtained 
from the exact estimation of cluster abundances (Liddle et al. 1995). The bound is weaker in 
the A-model than in the open model. 

4 Discussion 

A low-density universe is well motivated from several dynamical observations of galaxies and 
clusters. The simplest such low-density models are those in which the universe is open. In the 
context of inflation theory, however, we need a special idea such as the one-bubble inflationary 
scenario in order to produce an open universe. In this paper, motivated by the one-bubble 
inflationary universe scenario, we have examined the cosmological predictions based on the 
assumption that the scalar field is initially in the Bunch-Davies vacuum state. The initial 
perturbation spectrum has been derived by considering the evolution of perturbations in an 
open inflationary stage. Then the CMB anisotropics and the matter inhomogeneities have been 
examined. 

As the first test, we have performed a likelihood analysis for the CMB anisotropies by using 
the COBE DMR data. Interestingly, the COBE likelihood analysis gives severe constraints on 
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the model. Models with f^o < 0.4, 0,q < 0.5 are excluded at confidence levels of 92%, 83%, 
respectively. In a previous analysis associated with the conformal vacuum state (Gorski et al. 
1995), the likelihood function has another steep peak below f^o ^ 0.15. This complicates the 
statistical interpretation of the results (Gorski et al. 1995). In the case of the Bunch-Davies 
vacuum state, no such peak appears in the range of we are interested in, and so the likelihood 
analysis gives clear results. The COBE likelihood analysis is therefore a powerful probe of these 
open models. 

We have used the the COBE DMR maximum- likelihood normalization to predict the ampli- 
tude of matter fluctuations. According to this normalization method, there is little difference 
between the predictions of the Bunch-Davies vacuum and conformal vacuum cases. Even for 
the case Slo = 0.1, the discrepancy of as is a few percent. We obtain results that are similar to 
previous analyses: the power spectrum of the mass fluctuation fits the observations of galaxies 
and clusters for 0.3 ^ ^ 0.5. The required bias is unacceptably high for fio ~ 0.1, while high 
values of Qq demand anti-biasing. For example, f2o ^ 0.6 needs anti-biasing when h = 0.65. 
On the other hand, the A-models with Harrison Zel'dovich spectrum have higher amplitude 
compared with open models. The A- model needs anti-biasing for Q ^ 0.4 even when h = 0.65. 
The A- model therefore needs low h or a tilted spectrum (Ostriker & Steinhardt 1995). For the 
range 0.3 ^ Q 0.5, open models give unacceptably small bulk velocities compared with the 
POTENT analysis. However, given the present quality of the velocity data and the problem of 
cosmic variance, one might be reluctant to draw strong conclusions from this fact. The rough 
estimation associated with the galaxy formation gives lower bound of 0,q consistent with the 
value discussed above. 

It is very interesting that the COBE likelihood analysis has given the most severe constraint 
on this open model. The COBE likelihood analysis strongly prefers a high value of S7o- The 
peak value is around 0.7 55 0.8, and we can state that f^o > 0.5 with 83% confidence in 

this model. Considering both the COBE analysis and the matter inhomogeneity, we are led to 
prefer a value of Qq — 0.5 if the one-bubble inflationary scenario is correct. Such a model is 
consistent with the Press-Schechter analysis of the epoch of galaxy formation and is marginally 
consistent with the bulk velocity data. As the CMB data continue to improve, particularly on 
degree scales, we should be able to test this model. 

It is premature to rule out low inflationary models on the basis of this investigation 
at present, because we do not include the effect of bubble nucleation in the calculation of 
initial density power spectrum. Previous analysis indicates that the bubble nucleation effect 
in general excites fluctuations, and amplifies the perturbations on scales larger than curvature 
scale (Yamamoto, Tanaka, & Sasaki 1995; Hamazaki et al. 1995). One might therefore expect 
low-density models to fit the COBE data even more poorly once this effect is taken into account; 
however, since the calculation has not been done, we cannot be certain. In particular, the status 
of the super-curvature mode is still quite uncertain. 

Various modifications of the open model may also be viable. We must investigate the effect 
of gravity waves in an open inflationary universe. One might also consider the effect of tilting the 
primordial power spectrum; however, in order to improve the fit to the data one would probably 
need to tilt the power spectrum to have increased power on small scales, and such "blue" power 
spectra are not naturally produced by inflation. Such a model is probably too contrived to be 
plausible. 
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Table 1 

Confidence levels for Bunch-Davies open model 



fio(>) 


confidence level (%) 


0.1 


99.4 


0.2 


98.4 


0.3 


96.2 


0.4 


91.8 


0.5 


83.2 



Table 2(a) 

Amplitude of density perturbation for Bunch-Davies open model 



n 


A 






0"8 






(/i^Mpc) 3 




h=0.5 


h=0.65 


h=0.8 


0.05 


1.96(Q/27.9^K) 2 


x 10 2 


0.0041 


0.011 


0.021 


0.1 


2.42(Q/28.8^K) 2 


x 10 2 


0.032 


0.063 


0.099 


0.2 


2.83(Q/27.8^K) 2 


x 10 2 


0.15 


0.25 


0.35 


0.3 


2.94(Q/25.8//K) 2 


x 10 2 


0.31 


0.48 


0.66 


0.4 


2.90(Q/23.4^K) 2 


x 10 2 


0.49 


0.74 


1.00 


0.5 


2.72(Q/21.1^K) 2 
2.43(Q/19.3//K) 2 


x 10 2 


0.69 


1.01 


1.33 


0.6 


x 10 2 


0.89 


1.26 


1.63 


0.7 


2.03(Q/18.3AtK) 2 


x 10 2 


1.05 


1.45 


1.90 


0.8 


1.54(Q/18.3^K) 2 


x 10 2 


1.20 


1.66 


2.10 


0.9 


1.00(Q/18.9^K) 2 


x 10 2 


1.31 


1.80 


2.25 



15 



Table 2(b) 

Amplitude of density perturbation for conformal vacuum open model 



n 


A 

(/i^Mpc) 3 






h=0.5 


0-8 

h=0.65 


h=0.8 


0.05 


2.45(Q/18.7 / uK) 2 


X 


10 2 


0.0045 


0.012 


0.023 


0.1 


2.54(Q/23.1 / uK) 2 


X 


10 2 


0.032 


0.064 


0.10 


0.2 


2.89(Q/26.5/xK) 2 


X 


10 2 


0.15 


0.25 


0.36 


0.3 


3.01(Q/25.9^K) 2 


X 


10 2 


0.31 


0.49 


0.67 


0.4 


2.96(Q/23.5 / uK) 2 


X 


10 2 


0.50 


0.75 


1.01 


0.5 


2.77(Q/20.6 / uK) 2 


X 


10 2 


0.69 


1.02 


1.34 


0.6 


2.46(Q/18.3 / uK) 2 


X 


10 2 


0.88 


1.27 


1.64 


0.7 


2.04(Q/17.2 / uK) 2 


X 


10 2 


1.06 


1.50 


1.90 


0.8 


1.55(Q/17.5^K) 2 


X 


10 2 


1.20 


1.67 


2.10 


0.9 


i.oo(Q/i8.7//K) 2 


X 


10 2 


1.31 


1.80 


2.25 



Table 2(c) 

Amplitude of density perturbation for Harrison-Zel'dovich A-model 



^0 




0"8 






h=0.5 


h=0.65 


h=0.8 


0.1 


0.15 


0.29 


0.46 


0.2 


0.41 


0.68 


0.98 


0.3 


0.65 


1.0 


1.4 


0.4 


0.85 


1.3 


1.7 


0.6 


1.2 


1.7 


2.2 



Table 3 

Large scale bulk velocity for Bunch-Davies open model 



Q 




v R=40h -i Mpc (km/s) 






h=0.5 


h=0.65 


h=0.8 


0.05 


6.3 


8.6 


11 


0.1 


21 


29 


37 


0.2 


71 


90 


106 


0.3 


130 


160 


180 


0.4 


200 


230 


260 


0.5 


260 


300 


330 


0.6 


320 


360 


390 


0.7 


370 


410 


450 


0.8 


410 


450 


490 


0.9 


440 


480 


510 
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Figure Captions 



Figure 1. Power spectra of the CMB temperature anisotropy 1(1 + 1)C/ x 10 10 /27r for (a) the 
Bunch-Davics vacuum open model with f2o = 0.1, 0.2, 0.3, 0.4, 0.5. The data were provided by N. 
Sugiyama. These theoretical curves are normalized by the COBE likelihood. The curves are in 
descending order of Qq as one moves down at I = 50. The results of several CMB experiments are 
also shown, taken from the paper by Scott, Silk, & White (1995). To compare with degree-scale 
observations, careful investigations are required (Ratra et al. 1995). 

Figure 1(b). CMB power spectra for the conformal vacuum open model (Gorski et al. 1995). 
The curves and points are as in Figure 1(a). 

Figure 1(c). Power spectrum of the CMB temperature anisotropy 1(1 + 1)C/ x 10 10 /2-7r for the 
Harrison-Zel'dovich A-model with ^ = 0.1, 0.2, 0.3, 0.4, 1.0. (Sugiyama 1995) 

Figure 2. Contour plot of the likelihood function L(£Iq,Q) for the Bunch-Davies open model. 
The contour range is from L = 0.25 to L = 1.5, where the likelihoods are scaled so that L = 1 
corresponds to a flat Harrison-Zel'dovich spectrum with maximum-likelihood normalization. 

Figure 3. The marginal likelihood L m arg(^o) as a function of Qq f° r both the Bunch-Davies 
open model (solid line) and the conformal vacuum open model (dashed line). 

Figure 4(a). Power spectrum of density perturbation agP(k) for the open model with = 0.1, 
0.2, 0.3, 0.4, 0.5, for h = 0.75,0.70,0.65,0.65,0.60, respectively. We have taken n B h 2 = 0.0125. 
The points are from Peacock & Dodds (1994). 

Figure 4(b). Same figure as Fig. 4(a) but for the Harrison-Zel'dovich A-model with Qq = 0.1, 
0.2, 0.3, 0.4, 1.0. We have taken Q B h 2 =0.0128 and h = 0.8. 

Figure 5(a). Contours of <j(Mr = 10 12 M©, z = 1), in the (Qo — h) plane for the Bunch-Davies 
open models. The contour range is from a = 1.0 to a = 6.0. The dashed lines are a = 1.6 and 
2.4. 

Figure 5(b). Same figure as Fig.5(a) but for Harrison-Zel'dovich A-models. 
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